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Abstract: Exploring the current runoff characteristics after the large-scale implementation of the Grain 
for Green (GFG) project and investigating its sensitivities to potential drivers are crucial for water 
resource prediction and management. Based on the measured runoff data of 62 hydrological stations in 
the Weihe River Basin (WRB) from 2006 to 2018, we analyzed the temporal and spatial runoff 
characteristics in this study. Correlation analysis was used to investigate the relationships between different 
runoff indicators and climate-related factors. Additionally, an improved Budyko framework was applied to 
assess the sensitivities of annual runoff to precipitation, potential evaporation, and other factors. The 
results showed that the daily runoff flow duration curves (FDCs) of all selected hydrological stations fall 
in three narrow ranges, with the corresponding mean annual runoff spanning approximately 1.50 orders 
of magnitude, indicating that the runoff of different hydrological stations in the WRB varied greatly. The 
trend analysis of runoff under different exceedance frequencies showed that the runoff from the south 
bank of the Weihe River was more affluent and stable than that from the north bank. The runoff was 
unevenly distributed throughout the year, mainly in the flood season, accounting for more than 50.00% of 
the annual runoff. However, the trend of annual runoff change was not obvious in most areas. 
Correlation analysis showed that rare-frequency runoff events were more susceptible to climate factors. In 
this study, daily runoff under 10%—20% exceeding frequencies, consecutive maximum daily runoff, and 
low-runoff variability rate had strong correlations with precipitation, aridity index, and average runoff 
depth on rainy days. In comparison, daily runoff under 50%—99% exceeding frequencies, consecutive 
minimum daily runoff, and high-runoff variability rate had weak correlations with all selected impact 
factors. The sensitivity analysis results suggested that the sensitivity of annual runoff to precipitation was 
always higher than that to potential evaporation. The runoff about 87.10% of the selected hydrological 
stations were most sensitive to precipitation changes, and 12.90% were most sensitive to other factors. 
The spatial pattern of the sensitivity analysis indicated that in relatively humid southern areas, runoff was 
more sensitive to potential evaporation and other factors, and less sensitive to precipitation. 
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1 Introduction 


As an essential process of the hydrological cycle, the change in runoff plays a vital role in the 
sustainable development of human society, the economy, and the living environment (Greve et al., 
2014; Sun et al., 2016). In the past few decades, the runoff of 24% of the world's major rivers has 
experienced significant changes, most notably declining trends in Asia's large rivers (Li et al., 
2020). Yang et al. (2022) reported that from 1965 to 2018, the runoff of the Haihe River Basin, 
Liaohe River Basin, and Yellow River Basin in China showed a significant downward trend, 
leading to a sharp decline in the supply of ecosystem services and severe damage to human 
well-being. Therefore, it is of great significance to evaluate the influencing factors of runoff 
change for the development and utilization of water resources, agricultural production, and 
economic development (Sun et al., 2013). 

In order to improve the ecological environment, a series of ecological environment protection 
measures have been taken in the basin (Ouyang et al., 2013). Among them, the Grain for Green 
(GFG) project has dramatically optimized the land use structure and improved ecosystem services 
(Zeng et al., 2020). The study of Yang et al. (2022) showed that abrupt changes in runoff in many 
rivers occurred in the 1990s, indicating that the implementation of the GFG project can reduce 
runoff. However, climate change also drives the hydrological cycle, primarily through changes in 
precipitation and temperature (Xie et al., 2015). Under global warming, hydrological variables 
such as surface runoff, evapotranspiration, and precipitation change greatly (Nilawar and Waikar, 
2018). Therefore, studying the current runoff characteristics after the GFG project and 
investigating the runoff responses to climate-related factors will help us profoundly understand 
the water resource situation. 

The Weihe River Basin (WRB) is a typical transition watershed of arid and semi-arid areas in 
China. Its water resources are susceptible to climate change (Song et al., 2007). Yang et al. (2014) 
found that climate factors could promote runoff in Northwest China. Qiu et al. (2022) found that 
the magnitude and frequency of extreme precipitation in the WRB have strengthened in the past 
60 years and showed a strong aggregation mode. In addition, in some areas, higher potential 
evapotranspiration and temperature will also intensify the hydrological cycle (Li et al., 2016; 
Miao et al., 2016). As the largest primary tributary of the Yellow River Basin (YRB), the Weihe 
River not only provides a huge amount of water for the YRB and irrigates the Guanzhong Plain in 
the lower reaches of the Weihe River, but also maintains the ecological balance and social and 
economic development in Northwest China (Zuo et al., 2015). Since the large-scale 
implementation of the GFG project, runoff in the WRB has experienced a dramatic decreasing 
trend (Gao et al., 2013; Chang et al., 2015). 

The relationships between runoff and climate factors are usually assessed by hydrological 
models or empirical statistical methods in the WRB. Gao et al. (2013) used the double-mass 
curves for the period 1932-2008 to quantify the contribution of precipitation to runoff change, 
and the results showed that the impact of precipitation on runoff decline accounted for only 
17.20% before and after the transition year. Zhao et al. (2013) quantitively evaluated the impacts 
of climate change and human activities on runoff in the WRB from 1958 to 2008 based on 
Budyko curve method. They found that the decrease in annual runoff at most hydrological 
stations was due to the reduced precipitation and increased potential evaporation. Li et al. (2019) 
applied the Soil and Water Assessment Tool (SWAT) model in the WRB and found that in the 
tributaries and upstream hydrological stations, from 1970 to 2016, the dominant factor of runoff 
reduction was determined to be climate change, which accounted for more than 90.00%. 

Previous studies have focused on the long-term change in runoff and its driving factors in the 
WRB while ignoring the current runoff characteristics. In addition, there are significant spatial 
differences in the geographical characteristics and climate types of the whole WRB, which makes 
the response of hydrological regimes in different places to climate change quite different (Huang 
et al., 2016). Existing studies primarily analyzed the sensitivity of the runoff at the outlet of the 
watershed to the potential influencing factors. However, they ignored the detailed and systematic 
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analysis of the relationship between climate variables (i.e., precipitation, evaporation, and other 
related parameters) and runoff characteristics of the whole WRB. In this study, several classical 
hydrological analysis methods were applied to analyze the runoff characteristics and its 
sensitivities to climate factors in the WRB from 2006 to 2018. The aims of this study are to (1) 
investigate the temporal and spatial variations in runoff characteristics of 62 hydrological stations 
in the WRB; (2) estimate the correlation between selected runoff characteristics and climate 
factors; and (3) calculate the sensitivities of annual runoff to precipitation, potential evaporation, 
and other factors based on an improved Budyko method. 


2 Materials and methods 


2.1 Study area 


The WRB, located on the central Loess Plateau in northern China, is chosen as the study area. The 
main Weihe River has a total length of 818 km and a drainage area of 1.35x10* km? (Fig. 1). It 
originates north of Niaoshu Mountains in Gansu Province and eventually flows into the Yellow 
River at Tongguan County, Shaanxi Province. The topography of the basin decreases from west to 
east, with the variation range of elevation being 243-3916 m. There are many tributaries of the 
Weihe River. The Jinghe River and Beiluo River are the first and second-largest tributaries, 
respectively. Tributaries with a drainage area of more than 1.00x10? km? include Bangsha River, 
Jihe River, Heihe River, Fenghe River, and Bahe River in the south and the Xinhe River, Sandu 
River, Hulu River, Niutou River, Qianhe River, Qishui River, Shichuan River, Jinghe River, and 
Beiluo River in the north (Zhou and Yan, 2014). The WRB is located in the transition area 
between the arid and humid areas of China, with an average annual precipitation of 572 mm. It is 
characterized by the uneven distribution of more in the south and less in the north. The runoff is 
also characterized by uneven regional distribution, considerable interannual variation, and uneven 
annual distribution. The GFG project was first piloted in Gansu, Shaanxi, and Sichuan provinces 
of China in 1999 and was fully launched nationwide in 2002. Therefore, the large-scale 
implementation of GFG project in the WRB occurred after 2000. 
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Fig. 1 Topography, river networks, and hydrological stations of the Weihe River Basin (WRB) 


2.2 Data sources 


To ensure the time consistency and availability of runoff data measured by the selected 
hydrological station, we chose 2006-2018 as the study period. Daily runoff and precipitation data 
of hydrological and meteorological stations were acquired from Hydrological Year Book of the 
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Yellow River Basin (the Ministry of Water Resources of the People’s Republic of China, 
2006-2018). Other observational daily meteorological data, including daily maximum and 
minimum air temperatures of 22 national meteorological stations, were derived from the National 
Meteorological Information Centre, China Meteorological Administration 
(http://www.cma.gov.cn/). The daily precipitation and temperature data were added or averaged 
into monthly or annual values for subsequent calculation and analysis. 


2.3 Data processing 


2.3.1 Runoff characteristic indices 


Table 1 lists the 20 indices chosen for this study and their descriptions and calculations. These 
indices describe the characteristics of runoff from different aspects. Qı, Qs, Qio, Q20, Q50, Qso, 
Qoo, Qos, and Qoo refer to runoff at 1%, 5%, 10%, 20%, 50%, 80%, 90%, 95%, and 99% 
exceeding frequencies, respectively. Among them, some are related to high-frequency runoff 
events (Qo, Q90, Qos, and Q99) and some are related to low-frequency runoff events (Q1, Qs, Q10, 
and Q20). Extreme consecutive daily runoff indices include maximum daily runoff (Max1), 
maximum runoff for 7 consecutive days (Max7), maximum runoff for 30 consecutive days 
(Max30), minimum daily runoff (Mini), minimum runoff for 7 consecutive days (Min7), and 
minimum runoff for 30 consecutive days (Min30). Additionally, there are also indices showing low 
runoff variability rate (Raos-as0) and high runoff variability rate (Ras:as50). These runoff series can 
basically reflect the general hydrological characteristics of a watershed (Bassiouni and Oki, 
2013). 


Table 1 Description and calculation of runoff characteristic indices 


Index Unit Description and calculation 


Qi, Qs, Qio, Q20, Q50, Qso, mm/d Daily runoff under different exceeding frequencies. The subscript represents the 


Qoo, Qos, and Qog percentage of runoff exceedance probability. 

Min,, Minz, and Min39 mm Consecutive minimum daily runoff. The subscript indicates the consecutive days. 
Max,, Max7, and Max39 mm Consecutive maximum daily runoff. The subscript indicates the consecutive days. 
Ross:Q50 and Ras:950 / Low-runoff variability rate and high-runoff variability rate, respectively. Annual 


runoff is divided by the median annual runoff, providing an overview of how the 
flow duration curve changes over time. 


RC / Ratio of annual runoff to precipitation, reflecting the average water production 
capacity of the basin during a certain period. 


R mm Annual runoff. 


Note: Qi, Qs, Qio, Q20, Qs0, Qso, Qoo, Qos, and Qog refer to runoff at 1%, 5%, 10%, 20%, 50%, 80%, 90%, 95%, and 99% exceeding 
frequencies, respectively; Min, Min7, and Ming refer to minimum runoff for 1, 7, and 30 consecutive days, respectively; Max, Max7, 
and Maxo refer to maximum runoff for 1, 7, and 30 consecutive days, respectively; Rgos:a50 and Ras:aso refer to low runoff variability 
rate and high runoff variability rate, respectively; RC refers to runoff coefficient; R refers to annual runoff; "/" denotes dimensionless. 


2.3.2 Precipitation- and evaporation- related indicators 


To explore which climate factor is most closely related to runoff characteristics, we selected 
several indicators and their combinations to perform correlation analyses (Table 2). Among 
climate factors, precipitation is the most active climatic factor in the natural water cycle process, 
and many precipitation-related factors play an essential role in hydrological processes (Chang et 
al., 2015). Temperature mainly affects the hydrological process by affecting evaporation. Due to 
the strong seasonality of precipitation and evaporation, precipitation seasonality (Psi) and 
potential evaporation seasonality (Ep si) were calculated to express their seasonal characteristics. 
In addition, precipitation and evaporation usually affect runoff through joint action, so the 
combination of the two was considered. Among the eight selected indicators, three were related to 
precipitation (annual precipitation, Psi, and average rain depth), two were related to evaporation 
(annual potential evaporation and Ep si), two were related to the combination of precipitation and 
potential evaporation (aridity index and seasonal correlation between water supply and demand 
(CORR)), and one to temperature. 
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Table 2 Climate-related characteristics selected to estimate the runoff characteristics 


Predictor Unit Calculation 


ee 12 
Annual precipitation (P) mm P= XP : 
isl 


where i refers to the month; and P; is the monthly precipitation (mm). 
Precipitation seasonality (Psi) / Py is calculated as follows (Walsh and Lawler, 1981): 


P; =P") P-P/12|. 


Average rain depth (a) mm/d Mean precipitation during rainy days. 
Annual potential evaporation (Ep) mm __ &, is calculated by a modified version of the Hargreaves (Droogers and Allen, 
2002): 


E, =0.0013x0.408 x RA x (T, +17) x(TD-—0.0123P)"”* , 


where Tav is the average of the mean maximum temperature and minimum 
temperature for each month (°C); TD represents the difference between the 
average maximum temperature and minimum temperature per month (°C); and 
RA represents the extraterrestrial radiation (MJ/(m?-d)), which were acquired 
from Food and Agriculture Organization (FAO) of the United Nations. 


Potential evaporation seasonality / Ep si is calculated as follows (Walsh and Lawler, 1981): 


(Œ si - 
psi) E, =E, "DY [Epa -Ep /12 


> 


where E, and E,m are the annual and monthly potential evaporation (mm), 


respectively. 
Aridity index (¢) / j= E, 
P 
Seasonal correlation between / CORR is the correlation coefficient between monthly precipitation and annual 
water supply and demand (CORR) potential evaporation (Petersen et al., 2012). 


; 12 
Average annual temperature (TA) K TA= E +274.15)/12, 
i=l 


where T; is the monthly mean air temperature of the i" month (°C). 


Note: "/" denotes dimensionless. 


2.4 Methods 


2.4.1 Flow duration curve (FDC) method 


Flow duration curve (FDC) indicates the relationship between the magnitude and the frequency of 
daily, weekly, monthly, or yearly runoff for a specific watershed without considering the 
continuity of time, providing a time percentage greater than or equal to a certain flow occurrence 
in the whole time series (Cigizoglu and Bayazit, 2015). FDC is widely used in hydropower 
engineering design, water resource supply, irrigation planning, verification of hydrological model 
results in areas without data, and analysis of regional hydrological characteristics (Vogel and 
Fennessey, 1994). Based on long time-series daily flows, FDC is one of the most effective 
methods for analyzing the runoff distribution characteristics and variation in a basin (Gao et al., 
2015). FDC was applied in this study to reflect the variation characteristics of watershed runoff 
comprehensively and graphically during the whole study period. 

2.4.2 Correlation analysis 

Regression analysis was adopted to evaluate the correlation between runoff values and potential 
climatic factors listed in Table 2. The choice of linear, exponential, or power functions depends on 
the determination coefficient (R?). The fitting equation corresponding to the maximum R? is 
considered the optimal time series distribution. 

2.4.3 Trend test 

Mann-Kendall (MK) test is a nonparametric method (Mann, 1945; Kendall, 1990). Because this 
method does not require samples to obey a specific distribution, and the results are not easily 
disturbed by a few outliers, MK test is widely used in the trend test of hydrological series (Wen et al., 
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2019). Gao et al. (2017) introduced the principle and application of this method in detail. In this 
study, MK test was used to detect the trend of observed runoff series at different gauges in the WRB. 
2.4.4 Runoff sensitivity coefficient based on the Budyko framework 
Based on the Budyko framework (Budyko, 1974), Fu's function (Fu, 1981) stated that the average 
annual balance between evaporation and runoff could be presented as follows: 
E R 
F(¢,@) P 1 at (1) 

where E refers to actual evaporation (mm); P represents precipitation (mm); R represents runoff 
(mm); F is an analytical equation representing the evaporative fraction (E/P) or runoff ratio (R/P); 
@ is aridity index; and œ is a parameter that represents all other factors that affect the average 
annual distribution of precipitation (such as the type of soil and vegetation and topography). 

Rewriting Equation 1, whereby aridity index is replaced by potential evaporation/precipitation 
(E,/P) (Berghuijs et al., 2017): 


1 
R(P,Ep,@) =P Er, T E, o \o , (2) 
P P 


where Ep is potential evaporation (mm). 

The elasticity coefficient refers to the ratio of the growth rate of two interrelated indicators in a 
certain period (Zhang and He, 2016). Then the absolute elasticities of runoff to precipitation, 
potential evaporation, and other factors were obtained by deriving the partial differential 
expressions for each factor: 


1 
ORR (g+) 


RP T>? (3) 
OP/P +(+) 
ra] oO a 
ee ARR _ 9° (g +1) d 4) 
eTe g+(La gee 
1 oO a) 
wt mgg) n(o zi] 
POS OR/R _ g +1 o (5) 
Ro Galea ts 
-$+(14+4°)o 


where œ represents other factors; ER p, ég Zp, and ER are the absolute sensitivities of runoff to 
precipitation, potential evaporation, and other factors, respectively. 

To assess the relative importance of each influencing factor, the calculation of the relative 
sensitivity of runoff to these three factors is as follows: 


|En; 


6. = 2 (6) 

lee + lene, [+ leer 

where @, is the relative sensitivity of runoff to precipitation, potential evaporation, and other 

factors. 6, can vary from close to zero (i.e., almost unaffected by that specific factor) to close to 
one (i.e., almost completely affected by that specific factor). 


3 Results 


3.1 Temporal and spatial characteristics of runoff 


3.1.1 Exceedance frequency for daily runoff 
Figure 2 shows the flow duration curves (FDCs) of daily runoff for 62 hydrological stations in the 
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WRB during 2006-2018. The FDCs of the selected hydrological stations fell in three narrow 
bands with mean annual runoff values of 3.43—13.33 mm, 20.92—109.30 mm, and 229.03—507.25 
mm, represented by red, green, and blue lines, respectively. The results depicted significant 
variabilities in runoff in different areas of the WRB. Among them, the green curves accounted for 
the highest proportion, showing that the runoff variability at most stations was in the middle 
range. Between Q20 and Qgo, the slopes of FDCs were relatively gradual. The value of Qso varied 
from 0.00 to 0.79 mm. For high-frequency runoff events (Qoo, Qos, and Qo) and rare-frequency 
runoff events (Qı, Qs, and Qio), the FDCs of each station decreased as the percentage of 
exceedance increased. It can be seen that in the upstream area of each tributary, including the 
main stem of the Weihe River, Hulu River, and Malian River in the northern Jinghe River Basin, 
the runoff was low. The high-runoff areas were mainly located on the south bank of the WRB, 
including the Bahe River, Laohe River, Fenghe River, and Shitou River. 


100 


MAR (mm) (b) 


—3.43-13.33 
10 } —20.92-109.30 

A —229.03-507.25 

z 

: 
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Fig. 2 (a), flow duration curves (FDCs) of daily runoff in the WRB with three distinct ranges of mean annual 
runoff (MAR); (b), spatial distribution of hydrological stations with different ranges of MAR. Red, green, and 
blue sites are 6, 45, and 11 of the 62 stations, respectively, where MAR varied from 3.43 to 13.33 mm, from 20.92 
to 109.30 mm, and from 229.03 to 507.25 mm, respectively. 


3.1.2 Trend of daily runoff at different exceedance frequencies 


We divided the WRB into 62 subintervals based on the selected hydrological stations. The area 
controlled by a single hydrological station was divided into one interval; and the area between 
adjacent hydrological stations upstream and downstream of the river was divided into an interval. 
During the study period, the variation trends of Qos, Qso, and Qs in different intervals are shown 
in Figure 3. For high-frequency runoff events (Qos), 17 subintervals showed significant upward 
trends, which were located at the source of the WRB, the middle and lower reaches of the Jinghe 
River Basin, and the middle reaches of Beiluo River Basin; 2 subintervals showed a downward 
trend, and other subintervals showed no significant trend. For intermediate-frequency runoff 
events (Qso), nine subintervals showed a significant upward trend, most of them were located in 
the Jinghe River Basin; two subintervals showed a downward trend, and others showed no trends. 


4 Upward e@ Notrend Y Downward C] Basin boundary 


Fig. 3 Spatial distributions of runoff trends under high-frequency runoff events (Qos; a), intermediate-frequency 
runoff events (Q50; b), and rare-frequency runoff events (Qs; c) 
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For rare-frequency runoff events (Qs), three subintervals showed significant upward trends, and 
two showed a downward trend. In general, the runoff on the south bank of the Weihe River was 
more stable than that on the north bank, and high-frequency runoff events had more variability 
than intermediate- and rare- frequency runoff events. 


3.1.3 Seasonal distribution of annual precipitation and runoff 


Taking the observed runoff at Linjiacun, Xianyang, and Huaxian stations located in the upper, 
middle, and lower reaches of the main stream of the Weihe River, respectively, Zhangjiashan 
station in the Jinghe River Basin and Zhuangtou station in the Beiluo River Basin as examples, 
and taking the sum of observed runoff of Huaxian and Zhuangtou stations as the total runoff of 
the whole WRB, the monthly distribution of precipitation and runoff was analyzed (Table 3). For 
the selected hydrological stations, the monthly precipitation and runoff were distributed unevenly 
and mainly concentrated in the flood season (June—September) (Table 3). Precipitation in the 
flood season accounted for more than 60.00% of the average annual precipitation for the year. 
Runoff from June to September accounted for approximately 50.00% of the total flow volume. 
The results revealed the nonuniformity of the annual distribution of the WRB. 


Table 3 Distribution of monthly precipitation and runoff in the Weihe River basin (WRB) 


Proportion 
in flood 
season (%) 


Hydrological Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec 


Station element (mm) (mm) (mm) (mm) (mm) (mm) (mm) (mm) (mm) (mm) (mm) (mm) 


Linjiacun Precipitation 6.75 12.36 28.45 42.01 67.83 79.81 125.67 117.11 129.45 51.22 19.95 3.34 66.09 
Runoff 1.14 0.82 1.18 1.68 2.26 2.11 5.86 4.14 447 3.48 1.70 1.34 54.94 
Xianyang Precipitation 6.38 9.91 23.41 36.78 60.26 54.16 103.21 100.87 105.06 47.98 20.56 3.67 63.49 
Runoff 2.03 1.77 2.07 3.21 3.97 4.05 8.22 7.25 11.29 7.64 4.26 2.69 52.71 
Huaxian Precipitation 4.36 8.85 17.37 36.79 59.82 50.17 90.58 79.21 98.98 43.32 26.11 4.48 61.33 
Runoff 1.86 1.53 1.84 2.95 3.54 2.78 6.36 6.04 8.95 6.01 3.62 2.26 50.54 
Zhangjiashan Precipitation 6.48 9.48 23.07 34.05 47.83 56.15 88.88 75.76 95.82 38.58 18.12 3.34 63.63 
Runoff 0.97 1.11 160 1.28 1.18 1.21 405 3.83 3.33 2.29 1.68 1.23 52.27 
Zhuangtou Precipitation 4.20 8.80 13.24 29.63 50.93 54.05 87.05 77.34 82.67 42.42 20.36 2.19 63.68 
Runoff 0.82 0.99 1.68 1.39 1.16 1.06 3.19 2.99 3.12 2.18 1.51 1.16 48.75 
Total Precipitation 5.64 9.88 21.11 35.85 57.34 58.87 99.08 90.06 102.40 44.70 21.02 3.40 63.79 
Runoff 0.99 0.94 1.29 153 169 1.52 3.85 3.50 441 3.03 1.87 1.26 51.31 


3.1.4 Trend analysis of annual runoff in different subintervals 


The annual runoff trend analysis in the WRB is shown in Figure 4. The regions with significantly 
increasing inflow were distributed in the middle and upper reaches of the Beiluo River Basin, the 
upper reaches of the Jing River, and the south bank of the Weihe River. The regions with 
significantly decreasing inflow were distributed in the lower reaches of the Beiluo River and the 
upper reaches of the Weihe River. In most regions, the runoff showed no significant trend in the 
study period. 


3.2 Climate-related factors on runoff variability 


Taking the total runoff of the WRB, the effects of climate-related factors on runoff variability 
were analyzed. Table 4 lists the R? values between the runoff characteristics and precipitation- and 
evaporation- related control factors. According to the average value of R? sorted from high to low, 
the three principal factors affecting runoff were precipitation, aridity index, and runoff depth in 
rainy days, while the three least important influencing factors were precipitation seasonality, 
potential evaporation seasonality, and potential evaporation. In general, the factors related to 
precipitation have the greatest impact on runoff, while the factors related to evaporation are less 
informative. Each influencing factor performed differently toward every runoff characteristic 
index. The most important influencing factors for rare-frequency daily runoff events (Qi, Qs, Q10, 
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Fig. 4 Spatial distributions of annual runoff trends in different subintervals 


Table 4 R? values between runoff characteristics and climate-related factors 


Index P Py a E, Ep si @ CORR TA 
Qı 0.72" lin 0.22" ow 0.64" pow 0.20 pow 0.06 pow 0.69 pow 0.467 exp 0.29 pow 
Qs 0:72 bow 0.16* pow 0.70% tin 0.17 pow 0.14 tog 0.64 pow 0.347 tog 0.17 exp 
Qio 0.78* pow 0.14 * pow 0.63" pow 0.16 exp / 0.69 pow 0.347 iog 0.22 exp 
Qz 0.68 exp 0.06*pow 0.32" pow 0.25 exp 0.05*in 0.49 pow 0.20" og 0.42 exp 
Qso 0.06* pow 0.227 tin / 0.30 exp / 0.12 exp 0.067 tin 0.43 exp 
Qso / 0.184 tin 0.08 tin 0.06 tin / / 0.05 "tin 0.2 og 
Qoo / 0.15*tin / / / / 0.15 "tin 0.12 tin 
Qos / 0.07 "in / / 0.06 tin / 0.15*in 0.09*iin 
Qog / / / 0.07 ‘exp 0.07 tin / 0.19 "tin / 

Max; 0.58% exp 0.35* pow 0.48% exp 0.23 pow / 0.59 pow 0.50*in 0.37 exp 

Max; 0.81 exp 0.18" pow 0.66* pow 0.15 pow / 0.75 pow 0.38" in 0.23 exp 

Max3o 0.75* pow 0.18* pow 0.62* ow 0.17 pow / 0.72 pow 0.427 tin 0.26 exp 

Min, / / / / 0.15 tin / 0.22%iin / 

Min7 / / 0.04* pow 0.08* exp 0.07 tin / 0.11 iin / 

Min3o 0.07 tin / 0.13 tin / / / 0.05" tin 0.05 tin 

Ragos:as0 0.08* pow / 0.12" ow 0.06 pow / 0.19 * exp / 0.12* ow 

Ros:Q50 0.71 * pow 0.097 tin 0.75* pow 0.07 pow 0.12 tog 0.57 pow 0.28* pow 0.05 exp 
RC 0.404 tin 0.17 *exp 0.22* pow 0.25 pow / 0.27 jaw 0.267 tin 0.56 exp 

MAR 0.69* pow 0.13 * exp 0.45* pow 0.23-pow / 0.62 pow 0.35 "tin 0.36 exp 


Note: P, annual precipitation; Psi, precipitation seasonality; a, average rain depth; Ep, annual potential evaporation; E, si, potential 
evaporation seasonality; ¢, aridity index; CORR, seasonal correlation between water supply and demand; TA, average annual 


temperature. Different fitting functions were considered here. "Lin", "exp", "log", and "pow" represent linear, exponential, logarithm, 


and power functions, respectively. Only the values of R? with the strongest correlation are displayed. "/" denotes R?<0.05; "+" and "—" 
indicate positive and negative relationships between runoff characteristics and climate-related factors, respectively. 


and Q20), consecutive maximum daily runoff events (Maxi, Max7, and Max30), and high 
streamflow variability (Rqs:q50) were precipitation, aridity index, and runoff depth in rainy days. 
Influencing factors identified as important for runoff coefficient were precipitation and annual 
average temperature. However, no significant influencing factors were detected for intermediate- 
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and high-frequency runoff events (Qs0, Qso0, Q90, Qos, and Qos), minimum consecutive daily runoff 
(Min,, Minz, and Min30), and low streamflow variability rates (Ra95:050). 


3.3 Sensitivity analysis of runoff to its potential influencing factors 


3.3.1 Spatial hydroclimatic characteristics 

The spatial distribution characteristics of the average runoff ratio, aridity index, and other factors 
during 2006-2018 are shown in Figure 5. Different values were graded by natural discontinuities 
in ArcGIS 10.2. The runoff ratio in most areas of the WRB was 0.05-1.00, accounting for 
approximately half of the total hydrological stations. The value of aridity index was concentrated 
in 1.30-1.50, matching the climate conditions in the arid and semi-arid areas of the WRB. The 
spatial distribution characteristics of aridity index and runoff ratio were opposite. The runoff ratio 
was usually small in places with high aridity index, such as the northern WRB. Other factors 
showed the characteristics of high in the north and low in the south of the Weihe River. 
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Fig. 5 Distribution of spatial hydrological and climatic characteristics in the study area based on the Budyko 


framework. (a—c), spatial pattern of the runoff ratio, aridity index, and other factors, respectively; (d-f), values of 
the runoff ratio, aridity index, and other factors at different hydrological stations, respectively. 


3.3.2 Spatial pattern of absolute runoff elasticities 

Based on the obtained aridity index and other factors (Fig. 5), we depicted the spatial distribution 
of the sensitivities of runoff to precipitation, potential evaporation, and other factors in Figure 6. 
The stations with a sensitivity of runoff to precipitation (eg p) greater than 1.00 accounted for 
about 93.55% of the observed stations, suggesting that the change in runoff caused by 
precipitation was always equal to or larger than the change in precipitation itself. The average 
value of £g p was 3.85, indicating that an increase of 10.00% in precipitation would lead to a 
38.50% increase in runoff. Generally, relatively humid regions (aridity index<1.20) have lower 
Eg p Values (e.g., the southern region of the basin). The WRB is a typical arid and semi-arid 
transitional zone in China, and the runoff of this area has a far higher sensitivity to precipitation 
changes. In addition, ¢gp always exceeded the absolute sensitivity of runoff to potential 
evaporation (ER, Ep), Suggesting that the percentage change in runoff caused by precipitation was 
greater than that caused by potential evaporation. Compared with precipitation, potential 
evaporation has an opposite effect on runoff. The average £r was —2.85, indicating that an 
increase of 10.00% in precipitation would lead to a 28.50% decrease in runoff. Runoff in 
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relatively humid areas showed higher sensitivity to potential evaporation than in arid regions. The 
elasticity of runoff to other factors (ER œ) was similar to that to €R £p, with a mean value of —2.16. 
Li et al. (2021) proved that other factors was significantly correlated with climatic and artificial 
factors. €R Was generally lower in dryland regions (e.g., the north of the Beiluo River and Jinghe 
River Basin) and higher in relatively humid areas (e.g., the south of the WRB). The sensitivity 
analysis showed that precipitation could promote runoff, while potential evaporation and other 
factors inhibited runoff. 


(a) (b) (c) 
ERP N EREp ERw 
© 0.61-2.62 e —8.29- 4.30 e-3.26- -2.51 
© 2.62-5.15 e —4.30- -2.08 o-2.51- -1.78 
@ 5.15-9.29 e -2.08-0.39 e—1.78- -1.04 
C Basin boundary C Basin bo 


ERED 
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0 5 10152025 3035 40 45 50 55 60 65 0 5 1015 2025 30 35 40 45 50 55 60 65 0 5 10152025 30 35 40 45 50 55 6065 
Station number Station number Station number 


Fig. 6 Absolute sensitivities of runoff to precipitation (erp; a), annual potential evaporation (€R.Ep; b), and other 
factors (€R,«; c) in the WRB. (d-f), values of er,p, ER,Ep, and éR,o at different hydrological stations, respectively. 


3.3.3 Relative sensitivity of runoff to the precipitation, potential evaporation, and other factors 
The relative sensitivities of runoff to precipitation (6p), potential evaporation (@gp), and other 
factors (@) are showed in Figure 7. For 87.10% of the selected hydrological stations, precipitation 
was always the most critical contributor to runoff changes (6p>@gp and 6»), while the other factors 
were the dominant factor for the remaining 12.90% of the selected hydrological stations (@,>6p 
and Ogp). There was no land area that potential evaporation accounted for the most part. 


4 Discussion 


4.1 Impact of climate-related factors on runoff variability 


Runoff is highly variable in different intervals of the study area, which may be related to the local 
differences in soil texture, landform, runoff mechanism, watershed size, and river network density 
(Rossi et al., 2016). For a specific watershed, the landform and soil properties are relatively stable 
in a short period, and the vegetation cover will not change significantly. Therefore, meteorological 
elements have become essential factors affecting hydrological processes (Liu et al., 2016). Among 
them, precipitation plays a vital role in the hydrological cycle process. The results of the 
correlation analysis showed that precipitation-related factors greatly impacted runoff change (Table 
4). To show the relationship between daily runoff and daily precipitation more clearly, we 
constructed the FDCs for daily precipitation at 62 meteorological stations (Fig. 8). Because there 
are many days with no rainfall in a year, the diagram was drawn only for rainy days. There was 
substantial overlap among precipitation events when mean annual runoff varied across 3.00 orders 
of magnitude (Fig. 2), suggesting that precipitation among stations was relatively stable. 
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Fig.7 Relative sensitivities of runoff to precipitation (@), potential evaporation (@p), and other factors (@.) at 
different hydrological stations 
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Fig.8 FDCs for daily precipitation of 62 meteorological stations in the WRB 


The shape of the curves showed that precipitation has a firm consistency with the rare-frequency 
runoff events, indicating that the influence of extreme precipitation events on rare-frequency 
runoff was more potent than that on intermediate- and high- frequency runoff events. 
High-frequency runoff events, such as Qos and Qv, are usually regarded as the base flow (Beck et 
al., 2015), which is relatively stable and not easily affected. 

Daily precipitation can also be described by the proportion of wet to dry days representing rain 
frequency (f) and average precipitation depth on rainy days (a) (Rossi et al., 2016). Based on the 
resulting values, we transformed mean annual precipitation, a, and finto grid maps by the Kriging 
method in ArcGIS10.2 (Hevesi et al., 1992) (Fig. 9). It showed that the areas with high mean 
annual precipitation were sometimes related to large a (e.g., high a in the middle Beiluo River), 
high f (e.g., high f in the Hulu River), or both (e.g., south bank of the Weihe River), which 
illustrated that mean annual precipitation is affected by extreme weather conditions. The 
precipitation in the south and east of the WRB was more abundant than that in the north and west. 
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Combined with the spatial distribution of runoff characteristics, this may be the reason why the 
runoff was greater in the south than in the north, more in the east, and less in the west (Fig. 5). 
However, Rossi et al. (2016) found that rainfall variability played a secondary role in runoff 
variability. The most important driving force for daily runoff variation may be the nonlinear 
transformation from precipitation to runoff, rather than the daily rainfall statistics. Additionally, 
precipitation also influences runoff generation by affecting the soil water content, vegetative 
cover, evapotranspiration, etc. (Rodriguez-Iturbe, 2000). 

Figure 9 f shows the spatial distribution of other climate-related factors. In the higher altitude 
areas, the CORR was high, while potential evaporation and annual temperature was low. Potential 
evaporation was large in the northern WRB, but annual temperature was not the highest. From the 
results of the correlation analysis (Table 4), CORR, potential evaporation, and annual temperature 
have poor correlations with runoff. Potential evaporation and annual temperature have adverse 
effects on runoff, while CORR has positive effects. 

The temporal variations in precipitation and potential evaporation are usually considered the 
leading causes of runoff and evaporation changes over time (Greve et al., 2014). Evaporation is a 
term balancing input precipitation, which may strongly affect the conversion of precipitation into 
runoff at a daily scale to adjust runoff variability (Samain and Pauwels, 2013). Aridity index, 
calculated by precipitation and potential evaporation, determined the spatial distribution of 
precipitation to runoff and evaporation (Budyko, 1974; Greve et al., 2014). In general, various 
climate factors are interrelated and interact with each other, leading to the variation in runoff. 
With the advancement of global warming, the atmospheric water content will increase, which will 
make further efforts on the Earth's water cycle (Huntington et al., 2016). Therefore, the 
spatiotemporal distribution of runoff would be consequently affected (Zhang and Dong, 2013). 


(a) (b) (c) 


MAP (mm/a) À a (mm/d) 
= High: 837.63 High: 8.86 High: 0.42 
Low: 407.36 = Low: 4.59 Low: 0.25 


——) 
(d) (e) (f) 
CORR Ep (mm/a) TA (K) 
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E... 0.35 = Low: 873.48 E... 277.94 


Fig. 9 Spatial distribution of mean annual precipitation (MAP; a), mean precipitation in rainy days (a; b), the 
ratio of wet to dry days (f; c), the seasonal correlation between water supply and demand (CORR; d), mean annual 
potential evaporation (Ep; e), and average annual temperature (TA; f) in the WRB 


4.2 Explanations for temporal and spatial variabilities of runoff 


In the WRB, the difference in runoff on a time scale was not significant, but it showed distinct 
regional characteristics in space. Runoff from the southern and western rivers was more abundant 
and stable than that from the northern and eastern rivers (Figs. 2 and 3). The reasons for this may 
be explained as follows. 

First, there is a difference between geological and geomorphic conditions. The north bank of 
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the Weihe River is densely distributed with river networks, most of which mainly flow through 
the hilly and gully areas of the Loess Plateau and Guanzhong Plain. The role of loose loess soil 
and sparse vegetation in the distribution and regulation of precipitation on the Loess Plateau 
makes the annual runoff in this region smaller, but the interannual variation is larger (Bai et al., 
2020). The alluvial plain of the Weihe River is very flat, with a high degree of water conservancy, 
belonging to the low water yield area (Wu et al., 2022). The tributaries on the south bank 
originate from the northern foot of the Qinling Mountains, with short sources, rapid flow, steep 
rise and fall of the river, and good conditions for runoff production and confluence, such as the 
Qingjiang River, Heihe River, Laoyu River, Dayu River, and Luofu River. Coupled with small 
rainfall loss and abundant recharge, the south bank of the Weihe River was considered a high 
water-yield area in the WRB (Gao and Feng, 2019; Wu et al., 2022). The southern WRB covers 
an area of 1.48x10° km?, with a total of 3.84x10° m? of self-produced water resources, accounting 
for approximately half of the total water resources in the Guanzhong area, which provides water 
for life, production, and ecology (Gao and Feng, 2019). 

Second, the different intensities of human activities. The intensity of human activities on the 
south bank of the Weihe River is weak, and the degree of development and utilization of water 
resources is relatively lower. The nonagricultural population accounts for only approximately 
10% of the total population in the region, and the per capita regional Gross Domestic Product 
(GDP), local fiscal revenue, and urban and rural residents' income are lower than the average 
level of Shaanxi Province of China (Gao and Feng, 2019). For the north bank, the population of 
the Guanzhong Plain is about 2.0x107, the local population density is as high as 800/km?, and the 
population urbanization rate is about 60.00% (Jia and Yang, 2017). Sustained population growth, 
urban scale expansion, and industrial structure upgrading all have a large water demand (Gao et 
al., 2013; Zhan et al., 2014; Wang et al., 2019). Additionally, the Guanzhong Plain has a long 
history of irrigation, which is the principal form of water consumption. Irrigation water on the 
Guanzhong Plain accounted for 48.65% of the total water consumption in 2019 (Shaanxi Province 
Department of Water Resources, 2019). 

Since the 1970s, a large-scale GFG project has been implemented in the WRB, which has 
significantly changed the land use/land cover conditions and improved the regional ecological 
environment (Li et al., 2017; Jiang et al., 2021). By increasing the soil surface roughness, slowing 
the runoff speed, increasing infiltration, and intensifying evaporation, runoff was restrained (Yang 
et al., 2021). Compared with the south bank, the GFG project has been implemented more widely 
in the loess area on the north bank and has achieved more significant results (Zhang et al., 2018). 
Some researchers have proposed that excessive vegetation restoration may be a poor choice in 
arid and semi-arid areas because local precipitation conditions are not considered, which 
exacerbates the shortage of water resources to a certain extent (Jiang et al., 2019). In this study, 
the annual runoff entered a relatively stable low-flow period, with no significant trend from 2006 
to 2018. For seasonal distribution, runoff usually accounts for a large proportion in the flood 
season with abundant rainfall, which is consistent with most studies (Gao et al., 2013; Xie et al., 
2022). The runoff coefficient of the basin was relatively low during the study period, varying 
from 0.05 to 0.11, which is smaller than the existing study (Feng et al., 2016), suggesting that the 
ability of precipitation to transform into the runoff is weakened. Therefore, large-scale vegetation 
restoration programs need effective and sustainable management. 


5 Conclusions 


Based on the measured daily runoff and meteorological data of 62 gauges in the WRB from 2006 
to 2018, we analyzed the runoff characteristics and the correlations between the runoff 
characteristic indicators and climate-related factors in this study. The modified Budyko 
framework was applied to evaluate the sensitivities of runoff to precipitation, potential 
evaporation, and other factors. The main findings are as follows. 

(1) Under different exceedance frequencies, the runoff of different hydrological stations in the 
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WRB varied greatly. The FDCs of the selected 62 hydrological stations fall in three narrow bands, 
with the mean annual runoff spanning approximately 1.50 orders of magnitude. The runoff on the 
south bank of the Weihe River was more abundant and stable than that on the north bank, which is 
the primary water source of the WRB. 

(2) Precipitation-related factors have more influence on runoff characteristics than evaporation- 
related factors. In this study, rare-frequency runoff events (Qi, Qs, Qio, and Q20), consecutive 
maximum daily runoff (Max), Max7, and Max30), and the ratio of high-runoff variability (Ras-.a50) 
were strongly related to precipitation, aridity index, and average precipitation depth on rainy days, 
while high-frequency runoff events (Qs0, Qso, Qvo, Qos, and Qos), consecutive minimum daily 
runoff (Mini, Min7, and Min30), and the ratio of low-runoff variability (Raos:950) had little 
correlation with the influencing factors. 

(3) Compared with potential evaporation and other factors, runoff was more sensitive to 
precipitation. In this study, the runoff of 87.10% of the hydrological stations was the most 
sensitive to precipitation, while the remaining 12.90% was dominated by other factors. Aridity 
index is an important factor affecting runoff characteristics. In places with lower aridity index, 
such as the south bank of the WRB, runoff was more sensitive to potential evaporation and other 
factors than to precipitation. Potential evaporation and other factors had adverse effects, while 
precipitation had positive effects on runoff. In general, runoff change is comprehensively affected 
by various factors. 
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